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Some aspects of the fluid mechanics of Liquid penetrant testing are considered. Penetration is represented by 
surface tension driven flow into defects of small defect width to depth ratio. Defect width is chosen so that both 
gravitational and non-continuum effects may be ignored. Penetration time is found to follow a Rideal- Washburn 
relation, in which 



yR cos 



where / is time. / defect depth, (X the dynamic viscosity, y the surface tension. R the defect width and 6 the 
contact angle. The proportionality constant, however, is shown to be strongly dependent on defect geometry and 
penetrant application procedure. The effect of slight fluid elasticity is shown to he negligible. 
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Introduction ( D ) * n<i raH * °^ penetration is independent of the orienta- 

tion of the defect, and 
Liquid penetrants are frequently employed to locate micro- (e) the surface tension is a function of the substances 

scopic (surface-aceessible) defects in solid, non-porous, involved and temperature only. 

materials. The basic technique involves 1 1 , 2, 3] 

Restrictions a and h can be used to define a range of 

(a) cleaning the material to be tested characteristic defect dimensions. For Liquids composed of 

(b) applying the penetrant simple molecules, we may estimate molecular spacing [4] as 

(c) removal of excess penetrant 

(d) application of a developer (e.g., an absorbent coat- A = 3.5 X l()~ 1() m (1) 

ing) 

(e) visual inspection, and Happel and Brenner [5], report continuum flows, such as 

(f) post inspection cleaning. Poiseuille flow, persist for Knudsen numbers (ratio of molec- 

ular spacing to some characteristic flow dimension (say, R) 

Success of the process depends on the ability of surface a s high as 0.01. With consideration of eq (1) then, restriction 

tension forces to draw the liquid into surface-accessible a implies that 

defects. In this report the identification and quantification of 

the factors which affect the extent and rate of penetration ^ > 35 x 10 _2 m (2) 

(and the somewhat analogous problem of rate of development) 

are considered. Restriction b governs the largest defect size of interest. It 

requires that surface tension effects based on the minimum 

defect size be much larger than gravitational effects based on 

w , ., ,. , . nil a defect depth (/). Mathematically 

We consider some aspects 01 the continuum Hutu mechan- r 

ics of surface tension driven flows. We further restrict our 

y 

flow situations to cases in which — >> pgl^ (3) 

R 

(a) the minimum defect dimension is much larger than 

the intermoleeular distance of the liquid, where y is the surface tension, p is the density, and g is the 

_______ acceleration of gravity. For an air- water interface, interpret- 



2. Scope of the Investigation 



1 Numbers enclosed in square brackets identify references. ing ">>" as "100X", and with / = 100/?, this implies that 
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R < 27 fim 



(4) Equation 6 is of course simply Poiseuille flow with solution 



In what follows then we shall deal with defect sizes of roughly 
10 -2 fxm to 10 2 /xm. Gravitational and non-continuum effects 
are ignored. Note further that our interpretation of eq (3) 
requires 



R 

I 



« 1 



(5) 



Restriction c removes the possibility that the surface tension 
might vary in some manner with the flow. We shall indicate 
where this assumption, though employed, is probably inap- 
plicable and more detailed investigation warranted. 

3. Fundamentals 

The salient features of surface tension driven flow into 
defects may be illustrated by flow through a~< open ended 
capillary tube, 2 as shown in figure 1. It is a straightforward 
exercise to show, for steady axisymmetric flow at finite 
Reynolds number, 3 that in the limit as R/l — > 0, the 
equations of motion become to order 1, 



4 



d 2 u ±du\ 

dr 2 r dr) dz 



dp 



(6) 





Penetrant 



2R 



Figure 1. Capillary flow of penetrant in a small tube. 



2 That is, for R/l — * 0, we initially treat the closed downstream end of the defect as 
if it were at t — °°, and had little effect. A more realistic approach, which accounts for 
the possibility of entrapped air is taken up in a later section. 

3 The same result eould be reached by allowing/?? = UR/v —* 0, where U is mean 
velocity and v kinematic viscosity. However this is not generally as good an 
approximation as R/l —* for the range of parameters considered here, and in view of 
equation 5 it is unnecessary. 



GR 2 

4/JL 



(1 



r 2 /R 2 ) 



vhere 



G = 



dp 

dz 



(7) 



(8) 



In eqs (6), (7) and (8), u is the axial velocity, \L the dynamic 
viscosity, r the radial coordinate and dp/dz the applied 
pressure gradient. 

The volumetric flow rate can be found as 



Q = 7rGR 4 /8fi 
For a constant area tube Q may be written as 

Q = 7TR 2 dz/dt 

so that 

dz_ GR^ 

dt 8/x 



(9) 



(10) 



(11) 



Ap 



2y cos 
R 



(12) 



yhere 6 is the contact angle. Clearly then over a lengths, 



_ Ap _2y cos 8 



Rz 



Using eqs (11) and (13) 



dz y cos 6 
4^1 



dt 



R/z 



(13) 



(14) 



or integrating and using z = 0, at t = one finds for the 
penetration time, at z = /, 



t = 



2/y 

R cos Oy 



(15) 



Equation 15 is the Rideal- Washburn [6, 7] equation 4 and 
fundamentally identifies and relates the pertinent variables 
to the penetration time. Equation 15 has been experimentally 
confirmed by S. S. Kozlovskii [8] for liquid flow in glass 
tubes (with gravity) and within reasonable limits (10%) by 
Good and Lin [9]. Kozlovskii [8] notes correctly that equation 



As applied to horizontal capillaries or equivalently following from the choice of 
defect sizes in the present case for negligible gravitational effects. 
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For capillary flows, the surface tension is related to the 
pressure jump across the meniscus by [4] 



15, derived from a steady flow assumption, should only be 
employed for times 



t » pR 2 //JL 



(16) 



with 



provided/ satisfies 



v = f(x\* x 2 ) on r 



A 2 /' =2. 



with pR/jU being the viscous time scale For the range of R 
values considered here, however, eq (16) is not a serious 
limitation. Modifications of eq (15) to account for effects 
such as capillary narrowing and networks of capillaries are 
discussed in references [10-13]. If eq (15) is to be applied to 
a plug of solid powder (of interest in the development koff [15] and Langlois | 16]. 

process) a tortuosity factor (f) [14] must be introduced (i.e., We consider here the simple hut illustrative 

/ -* £/). The factor R/f; 2 may be evaluated experimentally. elliptical cross section. We adjust our coordinat* 

that the most general ellipse can be written as 



Detailed solutions for some specific cross secti< 
as additional mathematical details mav be found 



(22) 



(23) 

>ns, as well 
in Sokolni- 

case of an 
j system so 



4. The Effect of Non-Circular Cross Section 

An obvious extension of the theoretical results of the 
previous section is to a more realistic representation of 
defects. To this end we consider flow through non-circular 
cross sections. The problem of steady state flow through 
straight "pipes" of uniform but non-circular cross section is 
highly developed. Under certain simplifying assumptions the 
problem reduces to one in potential flow theory and further is 
analogous to the problem of torsion of an elastic bar | 15|. 
The development of the equations here closely follows 
Langlois [16]. 

Consider a flow, in the X3 direction of the form 



u :i — u(x u x 2 ); u u u 2 







(17) 



being driven at stead) state by a uniform pressure gradient. 

Because of the form of eq (17) the nonlinear, inertia terms 
are identically zero, and continuity identically satisfied, so 
that the momentum equation becomes 



A 2 u + G/fi = 0. 



Equation (18) is subject to 

u(x u x 2 ) 
on pipe periphery T. Here 



dx x 2 dx 2 2 



(18) 



(19) 



(20) 



Equation ( 18) is of course a Poisson equation for the velocity. 
The system of eqs ( 18-20) may be shown to be equivalent to 
the Dirichlet problem'' 



l Set 



A 2 n=0 



u(x u x 2 ) = G/2f4p{x u Kg) -/(*i, * 2 )1 



(21) 



+ 



*2 



= 1. 



The solution of eqs ( 18-20) mav be easily 
G a 2 b 2 



a 



2/u d 1 + Ir 
with volumetric flow rate 

ttG 



(24) 
found to be | L6] 
(1 - x 2 /a 2 ~ xl/b 2 ) (25) 



a 3 * 3 



4/u a 2 

Again G may be estimated as 

y cos 6 



P 



G = 



s/b 



(H 



(26) 



(27) 



and penetration time is found to be, at z = /. 



/ = 



2/ 2 /u. 



+ I,'" 



y cos {a+ b)ab 



(28) 



The ratio of penetration times for elliptic to circular pipes 
of equivalent area may be found from eqs ( 15) and (28) to be 



t e 



11. \Zab(a z + 6 2 ) 



t cir. 



ab(a + b) 



(29) 



Note that for b = 100a say, penetration into the ellipse takes 
9.9 times longer than into the equivalent area circle. 
Equation (29) is a reflection of the fact that a circular pipe 
produces more flow rate, at a given pressure gradient, than 
an elliptical one, at equivalent area. 

5. Entrapped Air 

It would seem possible, in applying penetrant, to trap air 
in some defects. The flow induced compression of this 
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trapped air would, by reducing the effective pressure gra- z/l — 1.5% (39) 

dient, reduce flow rate and increase penetration time. 

Let us consider the problem in which flow is developing so that the effect can be quite dramatic; i.e., the penetrant 

slowly due to a continuously changing axial pressure gra- will not penetrate 98.5% of the depth of the defect. 



dient. It may be shown that the condition 



(30) 



is sufficient to guarantee slow axial development. By main- 
taining finite Reynolds number, the equations of motion 
become 



d 2 u 1 du\ 
KM — + "- =\ P o 



dr 2 



r dr 



Pi + 



2y cos 6\ 1 



R 



I 



(31) 



Here P is the atmospheric pressure and P x the pressure of 
the trapped air (initially atmospheric). Let us further take 



vhere 



The volume V x is then 



that 



PiVt = c 



C = ttRHP 



ttR 2 (1 - z) 



Pi = 



Pol 



d-z) 
The velocity may be found to be 

[2y cos 6(1 - z) - PozR] 



with 



ii =/(r) 



f(r) 



IR(1 



4 M 



(1 



2 /R 2 ) 



(32) 



(33) 



(34) 



(35) 



(36) 



(37) 



Equations (36) and (37) clearly show that the flow develops 
slowly enough to retain the parabolic profile of Poiseuille 
flow while continuously adjusting its magnitude (and flow 
rate) in proportion to a constantly diminishing pressure 
gradient. 
From eq (36) we note that the flow will cease when 



z/l = 



2y cos 6 



RP + 2y cos 6 



(38) 



ForPo of 1 atmosphere (~1 X 10 6 dynes/cm 2 ), R of 10 2 
jLtm and an air-water interface, eq (38) shows that the flow 
will cease at 



The flow rate, here a function of z, can be found directly 



, x (2y cos 6(1 - z) - P zR) ttR 4 , x 

Q(z) = — ^ g =- Z -^— (40) 



IR(1 - z) 



8/i 



Penetration time may be found, by integration of eq (40) to 
be 



t = 



8/ulI 2 

PoR 2 



1 



2y_ 

PoR 



6 - 1 



i 2y a 

In cos 6 

P,R 



- JL - cos e - i - 
p R J i 

z/l 



2y_ 

PoR' 



cos 6 — 1 



(41) 



As one would anticipate, eq (41) shows that the time to 
reach a z/l given by eq (38) is infinite. Both eqs (39) and 
(41) indicate that trapped air can have a large effect on the 
penetrant process. 

It should be noted that unlike the case considered here, 
one major effect of the compression will probably be to alter 
the surface tension and/or change the shape of the interface. 
Additional research then, drawing on more detailed infor- 
mation on the physics of the interfacial region, may be 
fruitful. 

6. Effect of Fluid Elasticity 

This section is undertaken on the assumption that the 
addition of various dyes, fluorescents and emulsifiers might 
add a slight elasticity to the penetrating liquid. There is at 
present no experimental evidence to support (or refute) this 
hypothesis. However, as the effects of even slight elasticity 
on entrace type flows can be large [17], it would seem 
worthwhile to include at least the outline of an analysis. We 
should also note that perhaps "elasticizing" the fluid may 
prove to be useful in some aspects of the penetrant process 
(say in desensitizing the fluid to a certain defect size range). 

For simplicity, consider steady axisymmetric fully devel- 
oped pipe flow. For such a flow (termed viscometric), it is 
fair to state that the difference between a Newtonian and 
(slightly) viscoelastic fluid is that the latter will exhibit non- 
zero normal stresses. The creation of these extra normal 
stresses in entrance flows may lead to a large entry pressure 
drop [17, 18]. This extra pressure drop can substantially 
negate the applied pressure gradient. 
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We consider a fluid in which the shear and normal stresses P e -> P a (51) 

can be written as 

The effect of elasticity then is apt to be small. 
du 



c)r 



V- a ( 42 ) 7. Concluding Remarks 



= C a [ — ) (43) 



Some aspects of the continuum fluid mechanics of pene- 
trant flow have been examined. Characteristic defect sizes of 
dw\ 2 * 10~ 2 to 10 2 jLtm were considered. Gravitational and noncon- 

tinuum effects were neglected. The defect depth was every- 
where assumed much larger than crack size. 
We further take Fundamentally, the penetrant process was found to obey a 

Rideal-W ashburn [6, 7] type relation, so that penetration 
C a =C' a \fJL (44) time 

where A is a characteristic relaxation time. Equations (42- / 2 /x 

44) are not unrealistic for weakly viseoelastic fluids in pipe fty (>()S $ 

flow. 

Equation (42) taken together with more general equations Considerations of noncircular cross section and entrapped 
of motion [18] permit a Poiseuille solution. Then a j r- however, show thai the proportionality constant in 

equation 52 is substantially geometry dependent and possibly 

__ L-lR (45) a ' s<) dependent on the penetrant application procedure. 

dr 2{i dz Elastic effects were shown to be negligible. 

Geometric and process dependence, although somewhat 
expected, is certainly disquieting. Trapped air, for example, 
can essentially eliminate penetration ability oyer a certain 



that 



« = C' a \ — r 2 ( — I ) (46) range of defect sizes. Of course, some very fundamental 

V 9*/ / 

length, have not been considered. It would seem that a 
Using eq (46), it is algebraically straightforward [18] to c u .. . , . . . . 

& ' r ... caretully controlled experimental investigation, coupled per- 



A 

Vjl \czj / questions such as what constitutes a sufficient penetration 



find as an expression for the effective pressure gradient * .. , 

1 naps with computer simulations tor complicated geometries, 

R 

z z 8jli Vr 3 



n p \ u \ ma y he worthwhile. 

■Cl^-Pe'k). (47) 



Letting 
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